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Abstract 

The infrared divergences of QCD scattering amplitudes can be derived from an anoma- 
lous dimension F, which is a matrix in color space and depends on the momenta and 
masses of the external partons. It has recently been shown that in cases where there 
are at least two massive partons involved in the scattering process, starting at two-loop 
order T receives contributions involving color and momentum correlations between three 
(and more) partons. The three-parton correlations can be described by two universal 
functions Fi and f2- In this paper these functions are calculated at two- loop order in 
closed analytic form and their properties are studied in detail. Both functions are found 
to be suppressed like 0{m^/s'^) in the limit of small parton masses, in accordance with 
mass factorization theorems proposed in the literature. On the other hand, both func- 
tions are 0{1) and even diverge logarithmically near the threshold for pair production 
of two heavy particles. As an application, we calculate the infrared poles in the qq — > tt 
and gg tt scattering amplitudes at two-loop order. 
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1 Introduction 

In the past few years, much progress has been achieved in the understanding of the infrared 
(IR) singularities of massless scattering amplitudes in non-abelian gauge theories, which are 
characterized by an intricate interplay of soft and collinear dynamics. While factorization 
proofs guarantee the absence of IR divergences in inclusive observables [1], in many cases large 
Sudakov logarithms remain after this cancellation. A detailed control over the structure of IR 
poles in the virtual corrections to scattering amplitudes is a prerequisite for the resummation of 
these logarithms beyond the leading order [2-5]. Catani was the first to predict the singularities 
of two-loop scattering amplitudes apart from the 1/e pole term [6], whose general form was 
derived later in [7, 8] . An interesting alternative approach to the problem of IR singularities 
was developed in [9], where the authors exploited the factorization properties of scattering 
amplitudes along with IR evolution equations to recover Catani's result at two-loop order and 
relate the coefficient of the 1/e pole term to a soft anomalous-dimension matrix. Surprisingly, 
the explicit calculation of the two-loop anomalous- dimension matrix revealed that it has the 
same color and momentum structure as at one-loop order [7, 8]. 

In recent work [10], it was shown that the IR singularities of on-shell amphtudes in massless 
QCD are in one-to-one correspondence to the ultraviolet (UV) poles of operator matrix ele- 
ments in soft-coUinear effective theory (SCET) [11-13]. They can be subtracted by means of 
a multiplicative Z factor, whose structure is constrained by the renormalization group. It was 
argued that the simplicity of the corresponding anomalous-dimension matrix holds not only 
at one- and two-loop order, but may in fact be an exact result of perturbation theory. This 
possibility was raised independently in [14]. A non-trivial test of this prediction at three-loop 
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order was performed in [15]. Detailed theoretical arguments supporting this conjecture were 
presented in [16], where constraints derived from soft-collinear factorization, the non-abelian 
exponentiation theorem [17, 18], and the behavior of scattering amplitudes in two-parton 
coUinear limits [19] were employed to show that the anomalous- dimension matrix retains its 
simple form to three-loop order, with the possible exception of a single structure multiplying 
a function of two linearly independent conformal cross ratios built out of the momenta of four 
partons, which would have to vanish in all collincar limits. It was also proved that the coef- 
ficient of the Sudakov logarithm in the anomalous dimension obeys Casimir scaling at least 
through four-loop order. 

It is interesting and relevant for many physical applications to generalize these results to 
the case of massive partons. The IR singularities of one-loop scattering amplitudes containing 
massive partons were obtained some time ago in [20], but until very recently little was known 
about higher-loop results. In the limit where the parton masses are small compared with the 
typical momentum transfer among the partons, mass logarithms can be predicted based on 
coUinear factorization theorems [21, 22]. This allows one to obtain massive amplitudes from 
massless ones with a minimal amount of calculational effort. This method has been verified 
for two-loop QED amplitudes [22] and QCD amplitudes describing top-quark pair production 
at hadron colliders [23, 24]. A major step toward solving the problem of finding the IR 
divergences of generic two-loop scattering processes with both massive and massless partons, 
without the restriction to the limit of small masses, has been taken independently in [25] 
and [26] . Interestingly, one finds that the simplicity of the anomalous-dimension matrix does 
not persist at two-loop order in the massive case. Important constraints from soft-collincar 
factorization and two-parton coUinear limits are lost, and only the non-abelian exponentiation 
theorem restricts the color structures appearing in the anomalous-dimension matrix. At two- 
loop order, two different types of three-parton color and momentum correlations appear, whose 
effects can be parameterized in terms of two universal, process-independent functions Fi and 
/2 defined in equation (5) below. Apart from some symmetry properties, the precise form of 
these functions was left unspecified in [25, 26]. 

The goal of this paper is to calculate and study these functions at two-loop order. A brief 
account of our main findings was presented in [27] . In the following section we review known 
facts about the structure of the anomalous-dimension matrix governing the IR poles of scat- 
tering amplitudes in non-abelian gauge theories. Our main new contribution, the calculation 
of the universal functions Fi and /2, is described in Section 3, where we also analyze some 
properties of these functions such as their analytic structure and their behavior near threshold. 
Contrary to statements made in the literature, we find that Fi and /2 do not vanish in the 
limit where the velocities of two massive partons approach each other. We then discuss the 
particularly interesting limit in which the parton masses are small compared with the typical 
momentum transfers. We find that in this limit our results are compatible with factorization 
theorems proposed in [21, 22]. As an application of our general results, we present in Section 4 
the anomalous-dimension matrices and Z factors relevant for the parton scattering processes 
QQj99 ~^ + colorless particles, where we allow for the presence of e.g. electroweak gauge 
bosons or Higgs bosons in the final state. Using these results, we derive analytic formulae for 
all 1/e" pole terms at two-loop order in the squared qq —> tt and gg tt scattering amplitudes, 
some of which were so far unknown or only known in numerical form. We attach these results 
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in the form of a computer program. In Section 5 we study elastic quark-quark scattering in 
the forward hmit s,m^ ^ \t\, in which our general expression for the anomalous-dimension 
matrix reduces to the cross anomalous dimension studied a long time ago in [28]. This serves 
as a non-trivial check of our calculation. Our main findings are summarized in the concluding 
Section 6. 

2 General structure of the anomalous dimension 

We denote by |A1„(e, {p}, {m})) a UV-renormalized, on-shell n-parton scattering amplitude 
with IR singularities regularized in d = 4 — 2e dimensions. Here {p} = {pi, . . . ,pn} and 
{m} = (mi, . . . , nin} denote the momenta and masses of the external partons. The amplitude 
is a function of the Lorentz invariants Sij = 2aijPi ■ pj + iO and p? = mf, where the sign 
factor aij — +1 if the momenta pi and pj are both incoming or outgoing, and aij — — 1 
otherwise. For massive partons {rrii 7^ 0), we define 4- velocities Vi = Pi/rrii with vf = 1 and 
> 1. We further define the recoil variables Wij = —(lij Vi ■ Vj — iO. We use the color-space 
formalism of [29, 30], in which n-particle amplitudes are treated as n-dimensional vectors in 
color space. Tj is the color generator associated with the i-th parton and acts as an SU{N) 
matrix on its color indices. Specifically, one assigns {T^)af3 — f^/s for a final-state quark 
or an initial-state anti-quark, {T°')ap = —t^^a fo^ ^ final-state anti-quark or an initial-state 
quark, and {T°')i,c = —if"-^'^ for a gluon. We also use the notation Tj ■ T,- = T^T^ summed 
over a. Generators associated with different particles trivially commute, while = Ci is 
given in terms of the eigenvalue of the quadratic Casimir operator of the corresponding color 
representation, i.e., Cg — Cq — Cp — ^^^^ for quarks and Cg = Ca — N for gluons. Below, 
we will label massive partons with capital indices (/, J, . . . ) and massless ones with lower-case 
indices (i, j, . . . ). Note that color conservation implies the relation 

5]t, + 5]t, = o (1) 

i I 

when acting on color-singlet states. 

It was shown in [10, 16, 26] that the IR poles of such amphtudes can be removed by a 
multiplicative renormalization factor Z~^{e, {p}, {rn}-i A*)) which acts as a matrix on the color 
indices of the partons. More precisely, we have 

Z^\t, {p}, {zzi},/x) |M„(e, M, M))Lqcd^^^^ = finite (2) 

for e ^ 0. The quantity denotes the strong coupling constant in the effective theory, which 
is obtained after integrating out the heavy partons [26]. It is related to the coupling constant 
Q,QCD q£ QCD via the decoupling relation a^'^^ = ^a^. In the relation above, this matching 
relation must be applied to the scattering amplitude in full QCD, as indicated. To first order 
in as, the matching factor appropriate for rih heavy-quark fiavors in QCD reads [31] 

+ 0{a% (3) 
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with Tp — ^. Here mj denote the masses of the heavy quarks. Note that, as an alternative 
to (2), one can convert the expression for the Z factor from the effective to the full theory by 
replacing — a^^^. We will make use of this possibility in Section 4 to predict the IR 
poles of the qq — > ti and gg — > ti amplitudes in full QCD. 
The relation 

Z-\e, M, {m}, Pi) -f- Z(6, M, {m}, pt) = -T{{p}, {m}, pt) (4) 
- ampL — — 

links the renormalization factor to a universal anomalous-dimension matrix T, which governs 
the scale dependence of effective-theory operators built out of coUinear SCET fields for the 
massless partons and heavy-quark effective theory (HQET [32]) fields for the massive ones. For 
the case of massless partons, the anomalous dimension has been calculated at two-loop order 
in [7, 8] and was found to contain only two-parton color-dipole correlations. It has recently 
been conjectured that this result may hold to all orders of perturbation theory [10, 14, 16]. On 
the other hand, when massive partons are involved in the scattering process, then starting at 
two- loop order correlations involving more than two partons appear [25] , the reason being that 
constraints from soft-collinear factorization and two-parton coUinear limits, which protect the 
anomalous dimension in the massless case, no longer apply [26]. 

At two-loop order, the general structure of the anomalous-dimension matrix is [26] 

r(M, {m}, ^l)^Yl ^ ^-p(«^) 1^ + E 
(id) ^ 

- XI -^1^ 7cusp(/5/j, a,) + X l^i^s) + X T/ ■ Tj 7cusp(as) In 



+ }2 ir'^T^T'jT^KFMj.PjK.PKi) (5) 



{I, J) k '^IkVl Pk 

The one- and two-parton terms depicted in the first two lines start at one-loop order, while 
the three-parton terms in the last two lines start at 0{af). Starting at three-loop order also 
four-parton correlations would appear. The notation . . . ) etc. refers to unordered tuples 
of distinct parton indices. We have defined the cusp angles /J/j via 

cosh (3ij = ^^'^ = -(7: J vi -vj-iO^ wjj . (6) 

2m I rrij 

They are the hyperbolic angles formed by the time-like Wilson lines of two heavy partons. 
The physically allowed values for wjj are wu > 1 (one parton incoming and one outgoing), 
corresponding to > 0, or wjj < —I (both partons incoming or outgoing), corresponding 
to — —b + iir with real b > 0} The first possibility corresponds to space-like kinematics. 



^This choice implies that sinh/3 = — 1. Alternatively, we could have used Pu = b — in with 6 > 0, in 
which case sinh (3 = w \/\ — w~'^. We have confirmed that our results are the same in both cases. 
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in which case the universal functions in (5) are real. (For /2 phases can arise if some of the 
variables Vi ■ are time-like.) The second possibility corresponds to time-like kinematics, 
in which case these functions have non-trivial absorptive parts. For academic purposes it is 
sometimes useful to consider the analytic continuation of (6) to the Euclidean region — 1 < 
wij < 1, for which = ib with < 6 < tt is purely imaginary. The form factors are real 
in this region, which however does not correspond to a physical scattering process. Later 
in this work we will consider the large- recoil limit \wjj\ ^ 1, which is obtained whenever 
mjuij <^ In this limit we have = ln(2'u;/j) up to 0{l/wjj) power corrections. 

The physically allowed values of the kinematical variables wjj, wjk, and wki entering the 
function Fi are not arbitrary. To see this, note that in the rest frame of particle / (where 
1) we have 



vj -vk = {vi ■ Vj) {vi ■ vk) - cos5\/ {vi ■ VjY - 1\/ {vi ■ vkY - 1 , 



(7) 



where 5 is the angle between vj and Vk-, and all the scalar products satisfy Vi-Vj > 1 by defini- 
tion. Requiring that cos 5 be in the range [—1, 1] leads to the condition G{\wij\, \ wjk\-i \ wki\) > 
with 

G{x,y,z) = l + 2xyz - {x^ + y'^ + z^) . (8) 

Since in a three-parton configuration there is always at least one pair of partons either incoming 
or outgoing, at least one of the wu variables must be below —1, and hence the function Fi is 
expected to have a non-zero imaginary part. For /2 there is no non-trivial constraint on the 
allowed values of the kinematical variables, but also in this case there is at least one time-like 
invariant, so a non-trivial imaginary part arises. 

The anomalous-dimension coefficients 7cusp(cts) and 7*(q;s) (for i = q,g) in (5) have been 
determined to three-loop order in [16] by considering the cases of the massless quark and gluon 
form factors. Of particular importance for our discussion is the cusp anomalous dimension for 
hght-like Wilson loops, whose two-loop expression is [33] 



as fOis\'' 

7cusp(Q!s) = \- I ) 

TT \ TT / 



67 TT^X 5^ 
36"l2j^'^"9^^''' 



+ 0{al) , 



(9) 



where n; is the number of massless quark fiavors. Similarly, the coefficients 7^(05) for massive 
quarks and color-octet partons such as gluinos have been extracted at two-loop order in [26] 
by analyzing the anomalous dimension of heavy-light currents in SCET. In addition, the 

velocity-dependent function 7cusp(/3, in (5) was derived from the known two-loop anomalous 
dimension of a current composed of two heavy quarks moving at different velocity [33, 34] . A 
recent reanalysis of this anomalous dimension has led to the compact form 



7cusp(/?, Ois) = 7cusp(Qis) P coth/? 



Li3(e' 



-^^) + /3Li2(e-^'^)-C3+^/3+^ 



(10) 



coth (5 



U,{e-'^) - 2/3 ln(l - e"^^) - (1 + /3) - - | 



+ 0{al) 
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In the limit of large cusp angle one finds 

7cusp(/9, a.) = 7cusp(«.) /9 + C»(e-'^) = 7cusp(«.) ln(2w) + C»(^) . (11) 

Our goal in the present work is to further explore the structure of the three-parton cor- 
relation terms in (5), which are parameterized by two universal functions: Fi{x,y,z), which 
describes correlations involving three massive partons, is totally anti-symmetric in its argu- 
ments, while f2{x,y), which parameterizes correlations between a pair of massive partons and 
one massless parton, is an odd function of its second argument. The main technical challenge 
is to calculate these functions analytically at two-loop order, which involves the evaluation of 
a complicated, non-planar two-loop graph. This is the topic of the next section. 



3 Calculation of Fi and /2 

Our strategy will be to first consider the function Fi. To do so, we calculate the two-loop 
anomalous-dimension matrix of the soft Wilson-line operator Og = Sy^^ (0) Sy^ (0) Sy^ (0) with- 
out imposing color conservation. This is important, since for three partons in a color-singlet 
state, color conservation would imply that f'^^T^T^T^ = _johcj,aj,h ^j,c ^ j,c-^ ^ g. The 
operator Og consists of three time-hke Wilson lines 

Sy{x) ^ P exp (^-ig dtv ■ A^^ix + tv)T''^ (12) 

along the directions of the 4-velocities of three massive partons. The anomalous dimension of 
this operator is given by 



r5(M) = - Ti • T2 7cusp(/3i2, as) + permutations 



1=1 



(13) 



+ 6ir^^ri"r|r|Fi(/3i2,/323,/33i) +c(«D . 



The function Fi follows from the coefficient of the 1/e pole in the bare matrix element of the 
operator. We will then obtain /2 from a limiting procedure. 



3.1 Color-space formalism for Wilson coefficients and operators 

At this point a comment is in order concerning the color-space representation of on-shell 
scattering amphtudes in full QCD, and of Wilson coefficients and operators defined in the 
low-energy effective theory, which as mentioned earlier is a combination of SCET and HQET. 
In the effective theory, n-jet processes are described in terms of an effective Hamiltonian [16] 

\nf{{p}, {m})) = J2 On{{p}, {m}, /^) \Cn{{p}, {m}, /.)) , (14) 

where the short-distance Wilson coefficients are vectors in color space, whereas the operators 
On act as matrices on the color indices. The sum extends over all operators and coefficients 
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that are color conserving in the sense of relation (1). The sum also extends over operators with 
different Dirac structure, but this point is irrelevant to our discussion here. Taking the inner 
product of ITC^) with a given color state (c | produces the color-stripped amplitude {c\7i^), 
which is a c-number in color space. 

The renormalized operators On{lj) are related to the bare operators via a renormalization 
factor Z, so that 

Oniip}, {m}, /x) = 0^^(e, M, {m}) Z(e, {p}, {m}, /i) . (15) 

The definition of the anomalous- dimension matrix in (4) then implies the renormalization- 
group equation 

d 



dlnfi 



On({p}, {m}, ij) = -On{{p}, {m}, ij) rn({p}, {z^}, ij) ■ (16) 



Similarly, from the fact that the effective Hamiltonian in (14) is scale independent it follows 
that ^ 

|Cn(M, {zii}, ^)) = r„({p}, {m}, ^) |C„({p}, {m}, /x)) . (17) 

As explained in [10, 16], the parton scattering amplitudes in full QCD are given by the on- 
shell parton matrix elements of the effective Hamiltonian. Denoting these matrix elements by 
double brackets ((...)), and using that on-shell parton matrix elements of the bare operators 
Qbare scalclcss in the effective theory and are thus given by their tree-level expressions, we 
obtain 

\Mn{e, M, {m})) = ((^n^) Z{e, {p}, {m}, /.) |C„(M, {m}, /.)) . (18) 

The tree-level matrix elements arc given in terms of products of on-shell spinors and polar- 
ization vectors and act as unit matrices in color space, {{Ol^^)) — {{On^'^)) 1. It thus follows 
that 

\Mn{{p}, {m}, n)) = lini Z-i(e, {p}, {rn}, n) \Mn{e, {p}, {m})) 

— e— >0 — — , ^ 

_ (19) 

= Y.^^oT))\Cn{{p}Am}.^i)) 

is the finite, subtracted amplitude introduced in (2). This quantity obeys the same evolution 
equation (17) as the Wilson coefficients. 



3.2 Calculation of Fi 

Consider now the vacuum expectation value of the soft operator Og defined earlier. The 
renormahzation factor Zg for this operator is determined by the condition that {{O^^'^^)) Zg 
be UV finite, and its anomalous dimension Tg is then determined from (4). Here ((O^^^)) 
denotes the bare vacuum matrix element calculated in HQET. Since now we are interested 
in the UV poles of these matrix elements, it is necessary to regularize IR singularities other 
than in dimensional regularization (see below). At two-loop order, the above renormalization 
condition gives 

[((0^-))(°)zf + ((O^-)) (1)^(1) + ((Or^))(2)]^^ . = , (20) 
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where here and below the superscripts in parenthesis refer in an obvious way to the order in 
the expansion in powers of as/47r. The tree-level matrix element is {{Os))^^^ = 1. The equation 

(2) 

above thus expresses the two-loop renormalization factor Zs in terms of two contributions, 

Zi'^ = - [((Or^))(^) + ((0^-))(i)z(i)] ^^p^,^^ . (21) 

(2) 

The function Fi is derived from the pole terms in Zg with totally anti-symmetric color 
structure, so we can limit the discussion to Feynman graphs involving the color generators of 
all three partons. Diagrammatically, the first contribution on the right-hand side contains the 
UV poles of the planar and non-planar two- loop graphs shown in the first row in Figure 1. The 
second contribution corresponds to the UV poles of the one-loop diagrams with a counterterm 
insertion, as illustrated in the second row of the figure. In the calculation of the UV poles we 
regularize IR divergences by assigning residual external momenta li to the Wilson lines, with 
Ui = —Vi ■ li > 0. While the individual contributions depend on the Ui regulators, their sum 
does not. Also, for concreteness we perform the calculation with three outgoing Wilson lines 
in the fundamental representation. Afterwards we replace the color matrices arising from the 
Feynman rules by t° — >■ T° to convert to the color-space formalism. For an incoming line the 
color matrix would get transposed, and in addition one would pick up a minus sign since the 
velocity in the corresponding heavy-quark propagator is reversed. As a result, in this case the 
correspondence would be (— t'*)^ T", in accordance with the rules given in [29, 30]. 
We find that the 1/e pole terms in the sum of all diagrams can be written as 

{I,J,K) 

-^tr"^T,^T^T,^Fi'\(3,2,P23,f33i) + ... , (22) 
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where the dots represent finite terms and terms involving less than three different color gener- 
ators. The double pole is multiplied by a symmetric color structure and receives contributions 
from the two-loop planar and one-loop counterterm diagrams only. It is given by the square 
of the one-loop anomalous-dimension matrix, as required by renormalization-group invariance 
[16]. We emphasize that the two-loop coefficient F^"^^ not only receives a contribution from 
the non-planar, triple-gluon graph depicted in the first diagram in Figure 1, but also from the 
planar two-gluon diagrams and one-loop counterterm graphs. In [25] a formal argument is 
sketched which suggests that the sum of the two planar two-loop graphs shown in the first row 
in Figure 1 has a vanishing coefficient multiplying the anti-symmetric color structure due to a 
cancellation between two divergent integrals. We find by explicit calculation that such a can- 
cellation does not take place. Instead, the contribution to Fi from the planar and counterterm 
diagrams reads 

^(2)pW+CT ^ 4 ^ ^^^^^^^ ^^^^^^^ ^^^^^^^ 

X + 2/3,jln(l - e-^*0 - Li2(e-^'^") + \ 

D 

where the indices I,J,K run over permutations of (1,2,3). Moreover, as we shall see in 
Section 3.5, important cancellations take place when the planar and non-planar contributions 
are combined. 

We next consider the diagram with the triple gluon vertex, whose color factor if"'^'^ 
is totally anti-symmetric. Apart from this factor, the contribution of this diagram to Fi is 
given by the two-loop integral 

dfiki dfik2 d'^h 



^IJK Wij Vk ■ ki 



(24) 



kl kl kl {v-i ■ ki - oji) {v2 ■ k2 - 102) (vs ■ k^ - U3) 
= - y-^j eijK I{wij, wjK, Wki) + UV finite. 
The non-planar contribution to Fi is 

7-1(2) non— planar 4 \ , , /r>c\ 

-^^1 ^ - 2^eijKl{wij,wjK,WKi) . (25) 

I,J,K 

Introducing Feynman parameters and carrying out the loop momentum integrals, we obtain 

7(^12, W23, W31) = (^23 ^"31 + ^^12) J [dx] [dy] XiUs {xiX2 + X2X3 + XsXi)'^ (26) 

X [xiivl + yl + 2W23I/2I/3) + X2{yl + yl + 2w2,iyzyi) + X2,{yl + yl + 2wi2yiy2)\ ^ , 
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where [dx] = dxi dx2 dxs S{1 — Xi — X2 — X3), and all integrals run from to 1. This result can 
be recast into the five-fold MeUin-Barnes representation 



I{wi2,W23,W3i) = 2{W23W31 + Wu) 



Y\dzi 



i=l 



{2w2zf''-\2w3iY''-\2w^2 



2^2-1/ 



X 



r(i - 2;2i) r(i - 2^2) 



r(-2^3) r(-^4) r(^i + z^) r(zi + z^) v{z2 - z^) r(z3 + z^) 

(27) 



V{zi + Z2 + Z3 + Z4 + Z^) 

X r(2;i + Z2 + Z4) r(z2 + zs + Z4) r(z2 + Z4 + z^) r(i - Z2- Z4- z^) . 



Decomposing the wjj variables in terms of exponentials of cusp angles, wjj = cosh = 
(cK/j + ajj)/2 with ajj = e^^-' , we can convert the factors {2wijY^'^ i^^^ powers of au by 
introducing three more MeUin-Barnes parameters. By applying Barnes' Lemmas repeatedly, 
we can then reduce the representation (27) to a three-fold one: 



7(^12, W2Z, W31) = 2(W23 W31 + W12) 



{2TTif 



dzi dz2 dz^ a-^o^^ a^} ^^^o; 



l-2zi -l-2z2 



*31 



X T{-z^ - z^) r(l + z^- Z3) T{-z^ + Z3) r(l + zr + Z3) 

X V\-Z2 - Z3) r\l + Z2- Zs) r\-Z2 + Zs) F^l + Z2 + Z3) . 



(28) 



The remaining integrals can be performed by closing the contours and summing up the 
residues. The resulting expression for I is rather complicated, but the totally anti-symmetrized 
sum needed in (25) turns out to be amazingly simple: 



^(2)non-pIa.ax ^ _1 ^ ejjK^h^KI COth (3kI ■ 



(29) 



I,J,K 



In dealing with the Mellin-Barnes representations we have used the program package MB [36] 
and associated packages found on the MB Tools web page [37]. We have checked the answer 
for this diagram numerically using sector decomposition [38]. We have also checked that for 
Euclidean velocities our result for the triple-gluon diagram agrees numerically with a position- 
space based integral representation derived in [25]. Combining all contributions, we finally 
find 

(2), 



Fi {Pl2:/323:/33l) ^ 3 X] ^IJK Q^I^Ij) I^KI COth. (3kI, 

I,J,K 

where we have introduced the function 
g{p) = coth P 



P' + 2P ln(l - e-2^) - Li2(e-^'^) + ^ 

D 



TT 



(30) 



(31) 



The constant term — 7r^/6 has been added by hand, so that g{P) vanishes for /3 — > oo. Its 
effect cancels in the anti-symmetrized sum over terms in (30). 
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3.3 Derivation of fi 

While the three-parton contribution described by Fi is interesting on general grounds, there are 
not many processes of phenomenological importance in which three massive, colored partons 
are produced in a high-energy collision. For instance, searches for heavy, colored superpartners 
at the LHC will most likely focus on the pair production of squarks and gluinos. Hence, the 
three-parton term proportional to the function /2 in (5) is of greater practical importance. 
This function can be obtained from the result (30) by writing = —0-23^2 ■ Ps/^ria, Wai — 
— (T31 v\ ■ pz/mz and taking the limit 777.3 — > at fixed vi ■ ps. In that way, we obtain 

/2 (Pi2: In Z^Hi^^] ^ 3 lini Fi(/3i2, M ■ (32) 

\ — (T31 Vi ■ Ps/ m3-»0 

(2) 

Starting from the expression for given earlier, we immediately derive the two-loop coeffi- 
cient 

H /5i2, In = -4c/(/3i2) In , (33) 

V -(Ti3 Vi ■ P3 ) -(Ti3 Vx ■ P3 

where g{f3) has been defined above. We believe it is not an accident that the function /2 
is linear in its second argument, but that this feature persists to all orders of perturbation 
theory. The reason is that the logarithm 

-q'23 V2 ■ P3 ^ -2(723 ^2 " P3 _ -2ai3 Vl " PS 
-(Ti3 Vi-P3~ H /i 

is really the difference of two divergent coUincar logarithms, and in order for the scale de- 
pendence to cancel between terms depending on one of the two logarithms, the dependence 
should be single logarithmic. 

3.4 Properties of the universal functions 

We finish this section by collecting some useful properties of the three-parton correlation 
functions. We first note that, at least to two-loop order, we can rewrite the above relations in 
the suggestive form 

Fl{(3l2, P23, P3l) = IY1 ^^J^ ^ IcuspipKI, ^s) , 

(35) 

r fa 1 -0-23^2-P3\ / ^ f ^1 -^23^2-^3 

/2 Pi2: In ) = -— g{pi2) 7cusp(Qis) In , 

\ -ai3 vi-ps/ 47r -0-13 vi ■ p^ 

where 7cusp(/3, Q!s) and 7cusp(ci;s) are the cusp anomalous dimensions entering the two-parton 
terms in (5), and at one- loop order 

7cusp(/?, ots) = 7cusp(Qi5) r{l3) , with r(/?) = (3 coth/? , (36) 

where 7cusp(cKs) has been given in (9). Whether a factorization of the three-parton terms into 
a cusp anomalous dimension times a function of another cusp angle persists at higher orders 
of perturbation theory is an interesting open question. 
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-4 -2 2 4 -4 -2 2 4 

W W 

Figure 2: Graphical representation of the functions r{P{w)) and g{(3{w)). Real parts 
are shown as solid lines, while imaginary parts, existing for w < —1, are shown as 
dashed lines. The regions w >1 and w < —1 correspond to space- like and time- like 
scattering, respectively. The shaded region —l<w<l corresponds to Euclidean 
velocity vectors. 



It is useful to have expressions for the functions Fi and /2 in terms of the recoil variables 
wij = cosh/3/j. These follow from 



r(/3H) 



w 



— - In (^w + Vw"^ — 1^ , 



w 



Inlw + \/vS^ — \ ) In 



Li2 



w 



TT 



2 

In^ [w + - 1 j - ^ . 



(37) 



These functions are real for space- like w > 1, while for time-like w < —1 they have disconti- 
nuities given by 



-lmr{(3{w)) = 9{-w - 1) 

TT 

-Im g{p{w)) = 9{-w - 1) 

71 



W 



w 



= In [4(^2 _ 1)] + 2 In (I - w + Vw^ - l) 



(38) 



Here w is always defined with imaginary part — iO, see (6). The function /2 has further 
discontinuities in the logarithm given in (34), provided that one of the scalar products vi ■ 
or V2 -ps is time-like and the other one is space-like. In Figure 2 we show the real and imaginary 
parts of r{P{w)) and g{P{w)) as functions of w = cosh/3. 

It is also interesting to expand about the zero- recoil point {w = 1, (3 = 0) and the threshold 
point {w = —1, f3 = iir). For the functions r and g near zero recoil, we find (with w = 



12 



1 + /3V2 + /3V24 + ...,(3>0) 



(39) 



To expand the functions Fi and /2 themselves, we first note that at zero recoil /3i2 = and 
/523 =11^ — Psi- Then, from (39) and the relations 



r(i7r — (3) — r{(3) — it: coth (3 , 

g{in -p)^g{p)- 2in coth /3 ln(l - e^^^) - (tt^ + 2inP) (coth /3 - 1) , 
we find that for Wi2 — > 1 



(40) 



;im F,(A„/3.3,/330 ^2.2 



7^2A(/33l)+^7^S(/53l) 



lim /2(/3i2,ln— 

^2^0 \ — (7i3 ■Ui • P3 

where 

A{(3) = (coth /3 - 1) (/5 coth l3 - 1) , 



■P3 \ ^ 



(41) 



S(/3) =coth2/3 



Li2(e-^^)-^ + /3^ 



+ 2 coth /3 



l-/3-^-ln(l-e-^^) 



+ 2/3. 



(42) 



We will make use of these results when studying elastic quark-quark scattering in Section 5. 

Near threshold, the expansion of the functions r and g is (with w — — 1 — 6^/2 — 6"^/24 + . . . 
6 > 0, and p = -b + m) 



y + l-^+3+0(53), 

7r^ + 2z7rln(26) ^ ^ 57r^\ tt^ + 2z7rln(26) - 5i7r ^ ^ ^ 



(43) 



Based on the symmetry properties of Fi(Pjj, Pjk, Pki) and /2 (/3jj, In _''-'k'"J-Pk \ ^g^g ^.q^. 
eluded in [25, 26] that these functions vanish whenever two of the 4-velocities of the massive 
partons coincide. Indeed, this seems to be an obvious consequence of the fact that Fi is 
totally anti-symmetric in its arguments, while /2 is odd in its second argument. This rea- 
soning implicitly assumes that the limit of equal velocities is non-singular, but is invalidated 
by the presence of Coulomb singularities in r(/3) and g{(3) near threshold. In order to see 
this, consider the limit where two massive partons 1 and 2 are produced near threshold, i.e., 
Wi2 ^ — 1. In this case we can define the relative velocity Vu = V1—V2 in the frame where the 
two partons have equal and opposite velocities (for mi = 7712 this is the rest frame of the pair). 
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and — 4(t(;i2 + l)/(tyi2 — 1) — b"^ + 0{b^). To subleading order in the relative velocity, we 
thus have 

L-.^o = -^1±^|!« + (2 + ^) + . (44) 

It is then straightforward to derive the following limiting expressions valid near threshold (i.e., 
for W12 ^ -1): 

2 

lim F^{P,2,P23,P^l)^-<Jl3^\[n' + 2inH'2\v^2\)]r\P^,)-ing'{Ps^)} 

x(p.-^)+0(.l), (45) 

Vl^i2| i^^aiy 

hm = ^ [n^ + 2inH2\v,2\)] ^ ■ + 0{al) . 

/3i2^i7r V -CF^?.v^ ■ v?.^ 47r^ \v^o.\ \m 



In both cases the scalar products of the unit vectors are nothing but cos 6*, where d is the 
scattering angle formed by the momenta of particles 1 and 3. In the first relation the primes 
denote derivatives with respect to f3^i. The second relation is recovered from the first one using 
(32) and noting that (— cxia) r'(/33i) — >• 1 for |/?3i| — >• 00. The above results are anti-symmetric 
in the parton indices 1 and 2 as required (note that /3i3 = /323), but they do not vanish in the 
threshold limit. On the contrary, they diverge logarithmically in this limit. 

3.5 Limit of small parton masses 

A particularly interesting limit is that of large recoil, where all the scalar products wjj become 
large in magnitude. According to the definition wu = —sij/{2mimj). this limit corresponds 
to mirrij — > at fixed sij. In this case, we find for the non-planar contribution in (29) 

^(2) non-pW ^ A ^^"^ ^^"^ l^"^ (46) 
3 W22, W3I Wv2 

eijK ^^^^^ [2H'^wij) [H'^wki) + 1) - ln(2«;^,)] + (^(^) • 



I,J,K 



Similarly, the planar two-loop plus counterterm contributions in (23) can be expanded as 

^(2)planar+CT = ^ 1^ W23 W31 

^ 3 'U;23 W2,l W12 

+ Yl ^iJK ^"3^^'"'^ [2 - 2 H'^wij) ( ln(2«;^,) + 2) (47) 

I,J,K ^^-^ L 



+ \n{2wKi) + Y - 1 
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Note that the leading terms in these expressions are unsuppressed, but they cancel in the sum 
of the two contributions. We then obtain 



I,J,K 



2 ln{2wKi) 



\n'{2wij)-ln{2wjj) + — - 



+ 



(48) 



which can of course also be obtained directly from (30). We see that Fi vanishes in the limit 
\wij\ — > 00. Similarly, in the limit of large recoil we have 



In 



-tT23 V2 ■ P3 

-0-13 vi ■ Pa 



w 



12 



71 



\n\2wu) - H'^wu) + ^ - o 

D 2 



12' 

(49) 

It follows from these expressions that the three-parton correlation terms described by Fi 
and /2 vanish like (mrmj/sij)'^ in the small- mass limit. This finding is in accordance with 
a factorization theorem proposed in [21, 22], which states that massive amplitudes in the 
small-mass limit can be obtained from massless ones by a simple rcscaling prescription for the 
massive external legs. From (5), (11), and (48) we conclude that corrections to this simple 
factorization theorem are suppressed by two powers of mimj/sij. More specifically, we find 
that in the small-mass limit (see also the corresponding discussion for the two-parton terms 
in [26]) 



T{{p},{m},iJ,) = T{{p},ii) + Y, 



Ci 7cusp(as) In — A-f\as 
nil 



+ 0(!!«). (50) 



where r({p},/i) is the anomalous-dimension matrix in the massless case, and Aj^{as) = 
j'^{as) — 7^(q!s) (and similarly for hypothetical, heavy color-octet partons) is the difference of 
the single-parton anomalous dimensions belonging to massive and massless quarks. The extra 
terms give rise to one-particle Zj factors for each massive parton involved in the scattering 
process. As explained in [26], in the limit of small parton masses the product Yli removes 
the IR poles in the ratio of the massive to massless amplitudes (without including loops 
of heavy partons^). We emphasize that such a simple procedure will no longer work when 
0{m'jmj/s'jj) corrections are included, since in this case novel IR poles with non-trivial color 
and momentum correlations arise. 



4 IR singularities in top-quark pair production 

As a first application, we employ our formalism to calculate the two-loop anomalous-dimension 
matrices for top-quark pair production in the qq ti and gg ti channels and use them to 
construct the IR poles in the virtual corrections to these processes at two-loop order. These 
anomalous-dimension matrices form the basis for soft-gluon resummation at the next-to-next- 
to-lcading logarithmic (NNLL) order for generic kinematics, i.e., not necessarily restricted to 
the threshold region. 

■^One can account for heavy-parton loops by applying the inverse of the decoupling relation (3) to the 
combined Z factor, see Section 4.2. 



15 



4.1 Anomalous-dimension matrices 



The first step is to derive the exphcit form of the anomalous-dimension matrix (5) in a given 
color basis for the partonic amplitudes (see, e.g., [5, 39]). We adopt the s-channel singlet- 
octet basis, in which the tt pair is either in a color-singlet or color-octet state. For the 
quark- ant iquark annihilation process qiipi) + qk{P2) tiips) + tjipi), we thus choose the 
independent color structures as 



kl 



(51) 



For the gluon fusion process Qaipi) + 9biP2) Uipz) + tjip^), we use the basis 



Ci) = S'^' Sij , |C2) = i/""^ {t% , |C3) = d^'' {t% . 



abc /j.c\ 



mbc /j.c\ 



(52) 



Here a, b, i, j, k, I are color indices. We find that the anomalous-dimension matrix for the qq 
channel can be written in the form 



F - = 



Cf 7cusp(«s) la^ + Cp 7cusp(/534, as) + 27^(0!^) + 27^(0;^) 



N 

+ 2 



/ N 1 (-Sl3)(-S24) .r, x' 

7cusp(as) In 7 ^ 2 7cusp(P34, as) 



+ 7cusp(as) In 



y-s) 

(-gl3)(-g24) 
(-Sl4)(-S23) 





1 



(53) 



2N 

1 



+ 4^^(/534) 



Cf 
2 

-N 



^ 
^ 2 



+ C?(aD , 



where s = S12 is the square of the partonic center-of-mass energy. The term proportional to 
g{P34) stems from the three-parton contributions 



/2 /534,ln- 



-Sl4 



)+/2(/334, 



In 



-g24 
-S23 



^ 

" 2 

-AT 



(54) 



With the help of the second relation in (35) this can be recast into the product of g{l3u) times 
a conformal cross ratio [14] of four momentum invariants. Similarly, for the gg channel we 
obtain 



99 



iV7cusp(a.) In — + Cf7cusp(/?34, "J + 27^(a,) + 2-^^{as 
A* 



+ 



N 



( ^ , (-Sl3)(-S24) ,n 

7cusp(a^) In — I _2 7cusp(P34, 



-s) m 



/O o\ 



1 
\0 1/ 



(55) 



+ 7cusp(q;s) In 



(-gl3)(--S24) 
(-Si4)(-S23) 



/O 

1 



1 

2 


\ 


N 


Af2-4 


4 




iV 




4 



+ ^5(/334) 



/ f o\ 

-iV 
\ /_ 



+ 0{a 
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When deriving these results we did not impose momentum conservation. Therefore, they can 
also be used for processes involving additional colorless particles such as electroweak gauge 
bosons or Higgs bosons. 

The above anomalous- dimension matrices are the ingredients needed for soft-gluon resum- 
mation at NNLL order for general kinematics, a topic which we leave for future work. An 
interesting limit, which is often discussed in the literature, is the threshold hmit s — > 4mf . For 
this purpose it is convenient to define the quantity (3t = \/l — 4,ml/ which is related to the 
relative 3-velocity Vti between the top-quark pair in the center-of-mass frame by \vtt\ = 2j3t. 
The kinematical invariants can be expressed in terms of Pt and the scattering angle 6 between 
partons 1 and 3 in this frame according to 



P34 — iTT — In 



1+A 
1- A 



si3 = S24 = -7^ (1 - Acos6l) , si4 = S23 = -7; (1 + Acos6l) . (56) 



Using the expansions in (43), we find that in the threshold limit jJt ^ the two-loop 
anomalous-dimension matrices reduce to 



qq 



N 

+ T 





in 






/ in 


+ in — 







/cuspv 




1 



(57) 



+ ^ [n^ + 2inla{4:Pt)] cos0 



2n^ 



^ 
-TV 



and 



^ 99 



7cusp(Q!s) <N [In — - in ] - Cp 



in 



1 + 7^^L(A) + 27n«.) + 27^(«.) 



+ 



7cusp(a.) + iTT - 1 ) - 7cusp(A) 



+ ^ [7r2 + 2i7rln(4A)] cos^ 



-10 
\ / 



/o 0' 

1 
Vo 1, 



+ 0(A) + 0{al) , 



(58) 



where the two-loop expression for 7cusp(«s) has been given in (9), and 



Nai 



/cuspvr-iy 27^2 



m 



n 



;T7r 2-- -1 + C: 



6 



(59) 



arises from the threshold expansion of the two-loop coefficient of the recoil-dependent cusp 
anomalous dimension in (10). 
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We stress the important fact that, as a consequence of the Coulomb singularities present 
in the function (?(/?), the three-parton correlation terms governed by /2 do not vanish near 
threshold. Instead, they give rise to scattering-angle dependent, off-diagonal contributions in 

(57) and (58). These off-diagonal terms were not considered in two recent papers [40, 41], 
where threshold resummation for top-quark pair production was studied at NNLL order. 
Explicit results equivalent to the threshold-expanded anomalous-dimension matrices (57) and 

(58) were given in [41]. We agree with those expressions up to the terms originating from /2. 
We leave it to future work to explore if and how the results obtained by these authors need 
to be modified in light of our findings. 



4.2 IR poles in the qq — ^ tt and gg — > tt scattering amplitudes 

The results (53) and (55) for the anomalous-dimension matrices allow us to construct the IR 
divergences in the virtual corrections to the partonic amplitudes for top-quark pair production 
at two-loop order. The key relation is that the IR poles in the partonic scattering amplitudes 
can be absorbed into a multiplicative renormalization factor, as shown in (19). Expanding 
this relation in powers of cts /47r yields 



|_^(l),sing^ = Z«|Al(f)), 

where we have written the UV-renormalized QCD amplitudes as 



(60) 



\Mn) = 47ra, 



«.9 , . ,(!)> /a 



471 " ' V47r 



)'l-Mf) 



+ 



(61) 



Note that to predict the IR poles at two-loop order, one needs the UV-renormahzed one-loop 
amphtudes to 0{e), since these terms combine with the 1/e^ piece of the one-loop Z factor to 
give a 1/e pole term. 

In evaluating the relations (60), one must make sure that the UV-renormalized amplitudes 
and the Z factors are expressed in terms of one and the same coupling constant. As mentioned 
in Section 2, one possibility is to express all results in terms of ag defined with five active flavors, 
by applying the decoupling relation (3) to the UV-renormalized QCD amplitudes. After doing 
so, any dependence on the number Uh of heavy-quark flavors drops out of the results for the 
pole terms. Alternatively, by applying the inverse decoupling relation to the coupling in 
the Z factors, one can express the perturbative series in terms of af^^ deflned in full QCD 
with six active flavors {uh — 1 for the top quark, and rii — 5 for the remaining quarks). One 
then recovers the full Uh dependence in the singular parts of the QCD amplitudes. At the level 
of the renormalization factor, the conversion to six active flavors is accomplished by adding 
an additional piece to the expression given in [10, 16]. At two-loop order, we flnd the general 
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result 



a 



QCD / p/ 

/ J- r 



+ 



a 



An 



QCD 



4e2 2e 



^ = 1 + ^^^ — ( ^ 



47r 



(w.nAj, 3 \ 



+ §(ro-2/5o) + 3g, 



2rp 



E 

i=l 



fM ^ 1 ill ^ 1 



1 2 A^^ 



H m ^ 



(62) 



The coefficients are defined via the expansion 



n>0 



47r/ 



(63) 



and similarly for the quantity F' = — 2(7^ 7cusp(q;s), where C, = for the qq channel, and 
Ci = Ca for the gg channel. The /2 term enters the two- loop 1/e pole via Fi/e in (62). 
We emphasize that in the /^-function coefficient /^o = y Ca — | TpUi and in the two-loop 
anomalous-dimension coefficients Fi and F'^ in (62) the number n/ of active flavors only includes 
the massless quarks, not the massive ones. The Uh dependence of the full-theory Z factor is 
contained entirely in the terms shown in the third line. 

The result (60) is an exact prediction for the IR poles of the partonic amplitudes at 
two-loop order, which can be tested against explicit loop calculations. In practice, however, 
one is interested mainly in the real part of the interference terms {Ain^lAin^), since it is 
these which are needed to calculate partonic cross sections. For this reason, we give results 
for the interference terms rather than the amphtudes themselves. For the specific case of 
tt production, the full results for both the qq and gg channels are rather lengthy and are 
included as a computer program in the electronic version of this paper. In what follows, we 
will define the color decomposition used at two-loop order and describe to what extent we 
can compare our results with those available in the literature. As explained below, the three- 
parton correlations proportional to /2 do not appear in the interference of the Born level and 
two-loop amplitudes, neither in the qq channel nor in the gg channel. 

For the qq tt channel, the result for the interference term between the Born and two-loop 
amplitudes can be decomposed into color structures according to [42] 



2 Re _ 2{N^ - 1) (^N^A'^ + 5" + ^ C« + Nm Df + Nriu Dl 



(64) 



To compute the IR poles in the color coefficients above, we evaluate the general relation (60) 
using (62) for the renormalization factor and (53) for the anomalous-dimension matrix. In 

addition, we need the finite parts of the UV-renormalized one-loop QCD amplitude up to 
0{e), decomposed into the singlet-octet color basis. We have obtained these through direct 
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calculation, using some of the master integrals computed in [43]. After enforcing momentum 
conservation, the color coefficients are functions of the invariants s = S12, ti = S13 = S24, rrit, 
and /i. We have verified that our results for the IR poles agree with the numerical ones from 
[42], with the analytic results for some of the color coefficients given in [44, 45], and with 
the results in the small-mass limit from [23]. In our case, all pole coefficients are available 
in analytic form. Since the Born-level qq — > tt amplitude is proportional to the color-octet 
structure in (51) and the three-parton correlations proportional to /2 enter the anomalous- 
dimension matrix (53) only in the off-diagonal terms, the contributions from /2 in the squared 
matrix element first appears at three- loop order. This was noted independently in [41]. 

For the gg —>■ tt channel, we follow [24] and decompose the interference term between the 
Born and two-loop amplitudes into color structures as 

2 Re \M^''^)gg = {N^ - 1) (^N^A^' + n B^' + + ^ 

+ N'm Ef + N'nu El + m Ff + F,^ + Gf + ^ Gl (65) 

+ Nn^Hf + Nmn, Hf, + NnlH^ + | /f + ^ /f. + f 4^) ■ 

The IR poles in the color coefficients are obtained as for the qq channel, except in this case we 
use the anomalous-dimension matrix (55). Results in the literature are available only in the 
small-mass limit [24] , and we have checked the agreement of our exact results with this limiting 
case. Since the exact results are new, we list in Table 1 the numerical values for the poles of 
the color coefficients at the point ti = —0.45s, s = 5mf, and ^ = m^. Again in this case the 
results do not depend on /2, the reason being that the corresponding contribution is multiplied 
by a color structure which is anti-symmetric under the exchange of the two initial-state gluons, 
while the gg — > tt amplitude is symmetric under this exchange. 



5 Elastic quark-quark scattering in the forward limit 

Another interesting application of our general formalism is the case of elastic quark-quark scat- 
tering at high energy and fixed momentum transfer (s, ^ |t|). The anomalous-dimension 
matrix for this case was analyzed at two-loop order in [28] by studying the cross singularities of 
self-intersecting Wilson loops. We will now show that the results derived in that paper can be 
obtained by taking a certain limit of our general results, and that this provides a cross-check 
on our calculation of the three-parton correlations governed by the function Fi. 

Consider the clastic process qij{pi) + q2iip2) —>■ Quips) + (l2k{pi) for massive quarks (mi = 
m2 = m) in the forward limit 

s, > > AjcD ■ (66) 
Here i,j,k,l are color indices, and 1,2 label the quark flavors. The relevant cusp angles can 
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(67) 



Table 1: Numerical results for the IR poles in the color coefficients (65) for top- 
quark pair production in the gg — > tt channel, evaluated at the point ti — —0.45s, 
s = 5m^, and n = rrit. The blank entries are not present in general, while the entries 
with value 0.0 vanish only for the particular choice /i = m^. 

be expressed in terms of the invariants s = (pi + ^2)^ and t — {pi — pa)^ as 

(S — 2?7T.^ \ 
2mP~) = ~ ^' 

Pi3 = P24 = arccosh(^^— 2— j = Oy—^j , 

- /523 = arccosh(^±^^) =^ + o{-^), 

where cosh 7 = Vi ■ V2 — pi •p2/m^- In the limit t/m^ — > these angles are described in terms 
of a single variable 7 > 0. Starting from the general expression (5), we then obtain for the 
cross anomalous-dimension matrix 

rcross(7) '^s) =^qq{s,t,1Tl )A')|_(^5 „j2 

• T2 7cusp(i7r - 7, a,) +Ti-Ts 7e„,p(0, a,) + Ti • T4 7cusp(7, «.)] (68) 
+ 47«(a,) + 24zr'' T^T'^T^ F,{m - 7, 7, 0) + 0{al) , 

21 



= -2 



where we have used color conservation to simphfy the various terms. The velocity-dependent 

cusp anomalous dimension 7cusp(/^! (^s) has been given in (10). Moreover, noting that the case 
at hand corresponds to the zero-recoil limit discussed in Section 3.4, we can read off Fi from 
(41): 



Fi(i7r -7,7,0) 



127r2 



7r^^(7) + mB{'^) 



(69) 



where A and B were defined in (42). 

To give explicit results for the anomalous-dimension matrix (68), we must first specify a 
color basis. We shall use the t-channel singlet-octet basis, where the two color structures are 



\ci)^5ij5ki, \c2)^t^jtli. 
In this basis, the result for the cross anomalous dimension valid to two-loop order is 
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r'cross(7) 
rcross (7) 
Tcross (7) 

cross V / / 



11 



12 



21 



0, 

N 



J 22 



N 



7cusp(7, ««) - 7cusp(«7r - 7, as) + dCpFi {in - 7, 7, 0) , 
7cusp(7, Ois) - 7cusp(«7r - 7, as) - 12iVFi(i7r - 7, 7, 0) , 
7cusp(«7r - 7, as) - 7cusp(7, "s)] + N 7cusp(7, "s) - 7cusp(0, ag) 



(71) 



We have simplified the diagonal matrix elements using that 4:'j'^{as) = — 2Cf7cusp(0, as), which 
follows from the expression for the anomalous dimension of a heavy-quark current derived in 
[26]. Expanding the cross anomalous dimension as (this conforms with the notation used 
in [28]) 

rcross(7, = ^ rw ,(7) + ( v)' r[-s(7) + ■ ■ ■ , (72) 

the explicit one-loop result is 



r(^) (7) 

cross V / / 







in coth 7 



2in coth 7 A^(7 coth 7 - 1) - ^ coth 7 ^ 



(73) 



whereas at two-loop order we obtain 
= 0, 



r(2) (7) 

cross V / / 

r(2) (7) 

cross V / / 

r(2) (7) 

cross V / / 



11 



12 



21 



n'^A{'y) + in ( ^(7) + 
31 lOTp 



31 5Tf 



18 9A^ 



36 9A^ 
rii 1 in coth 7 , 



ni 



coth 7 



22 



r(') (7) 

cross V / / 



22 



^' J ,2 

— { coth 7 



^2 ^3 



Li3(e-2^) + 7 Li2(e-2^) - (3 + y 7 + y 



(74) 



+ coth 7 



TT 



67 TT^ 



18 



Li2(e-2^) - 27ln(l - e'^^) - :-+(^-:-)7-72-^ 



TT^ 49 



5 



+ r + ^- T^^-n ^^i^^' (7C0th7 - 1) 



18 



9 



7rM(7) - «7r 5(7) + 



coth 7 I . 



"31 _ lOTp 

18 ~ "a/v" 

In deriving these expressions we have used the remarkable relation 
7cusp(7, Ois) - 7cusp(«7r - 7, ocs) = i7rcoth7 











TT 





36 



+ 6ArFi(i7r-7,7,0) + C>(a^), 



5 
9 



(75) 



which links the complicated 7-dependent terms in the difference of cusp anomalous dimensions 
to those in Fi{m — 7, 7, 0). In order to compare with the expressions given in [28], we must 
convert our results to the color basis 



kj 



Results in that basis can be obtained from ours by the rotation 

•'^ cross (7) '^s) — ^ rcross(7) '^s) ^ i 

where 



1 

1 



V-1 



1 i 

2 



(76) 

(77) 
(78) 



After performing this conversion, we find complete agreement with the results in [28].^ Since 
at two-loop order the off-diagonal elements in (71) receive contributions from Fi, this is a 
non-trivial check on our general result for this function. 



6 Conclusions 

The IR divergences of scattering amplitudes in non-abelian gauge theories can be absorbed into 
a multiplicative renormalization factor, whose form is determined by an anomalous- dimension 
matrix in color space. At two-loop order this anomalous-dimension matrix contains pieces 
related to color and momentum correlations between three partons, as long as at least two 

^The cross anomalous dimension calculated here is the transpose of that in [28] , due to a different ordering 
of the matrices in the renormalization-group equation (16) compared to Eq. (2.9) in that paper. Also, these 
authors set n; = 0. 
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of them are massive. This information is encoded in two universal functions: Fi, describing 
correlations between three massive partons, and /2, describing correlations between two mas- 
sive and one massless parton. In this paper we have calculated these functions at two-loop 
order. For Fi, this involved extracting the UV divergences of the vacuum matrix element of an 
HQET operator built out of three soft Wilson hues. The most comphcated technical aspect of 
the calculation was the evaluation of a non-planar two-loop diagram, which we accomplished 
using Mellin-Barnes representations. The function /2 was then obtained from Fi by taking 
the limit where one of the three partons becomes massless. 

Using the exact analytic expressions, we studied the properties of the three-parton correla- 
tions in the small-mass, zero-recoil, and threshold limits. We found that the functions Fi and 
/2 vanish as {mimj / siif in the small-mass limit, in accordance with existing factorization 
theorems for massive scattering amplitudes [21, 22]. On the other hand, and contrary to naive 
expectations, the two functions do not vanish in the threshold limit, where the velocities of 
two heavy partons become nearly equal. The reason is that Coulomb singularities arise in 
this limit, which compensate a zero resulting from the anti-symmetry under exchange of two 
velocity vectors. 

Our results allow for the calculation of the IR poles of an arbitrary on-shell, n-particle 
scattering amplitude at two-loop order, where any number of the n partons can be massive. 
As an application, we have derived the anomalous-dimension matrices for top-quark pair 
production (in association with colorless particles such as electroweak or Higgs bosons) in the 
qq — > tt and gg — > tt channels. These matrices form the basis for soft-gluon resummation at 
NNLL order for general kinematics, and in particular near the production threshold. We will 
explore in future work to what extent the new off-diagonal entries in the anomalous- dimension 
matrices, arising from the three-parton correlation terms, affect the numerical results for ti 
production at the Tevatron and LHC. Finally, we have used these matrices to determine the 
IR poles in the virtual corrections to the tt production cross sections at two-loop order in 
closed analytic form. The corresponding expressions, which agree with those in the literature 
where they exist, are very lengthy and are provided in the form of a computer program. 

As a second interesting application, we have studied the elastic quark-quark scattering 
amplitude in the forward limit, where s, ^ \t\. In this special case, our general expression 
for the anomalous-dimension matrix reduces to the cross anomalous dimension studied a long 
time ago in [28]. Since the three-parton correlation terms in the anomalous dimension give a 
non-zero contribution in this example, the fact that we find full agreement with the two-loop 
expressions given in that paper provides a non-trivial check of our calculation. 

The present paper completes the study of IR divergences of two-loop scattering amplitudes 
with an arbitrary number of massive and massless external particles, and in arbitrary non- 
abelian (or abelian) gauge theories with massless gauge bosons. In spontaneously broken gauge 
theories, our results can still be used above the symmetry-breaking scale, where gauge-boson 
masses can be neglected at leading power. At the symmetry-breaking scale a matching is done 
onto a non-interacting theory, in which the massive gauge bosons are integrated out. In this 
way, it is possible to resum electroweak Sudakov logarithms using effective-theory methods, 
as worked out in detail in [46] . The next step should now be to apply these general results to 
specific hadron-coUider processes. 
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